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Abstract 

The flux trapping that occurs in the process of cooling 
down of the superconducting cavity is studied. The critical 
fields B c 2 and B c 1 depend on a position when a material 
temperature is not uniform. In a region with T cx T c , B c 2 
and B c i are strongly suppressed and can be smaller than 
the ambient magnetic field, B a . A region with B c 2 < B a 
is normal conducting, that with B c \ < B a < B c 2 is in 
the vortex state, and that with B c \ > B a is in the Meissner 
state. As a material is cooled down, these three domains in¬ 
cluding the vortex state domain sweep and pass through the 
material. In this process, vortices contained in the vortex 
state domain are trapped by pinning centers distributing in 
the material. A number of trapped fluxes can be evaluated 
by using the analogy with the beam-target collision event, 
where beams and a target correspond to pinning centers and 
the vortex state domain, respectively. We find a number of 
trapped fluxes and thus the residual resistance are propor¬ 
tional to the ambient magnetic field and the inverse of the 
temperature gradient. The obtained formula for the resid¬ 
ual resistance is consistent with experimental results. The 
present model focuses on what happens at the phase tran¬ 
sition fronts during a cooling down, reveals why and how 
the residual resistance depends on the temperature gradi¬ 
ent, and naturally explains how the fast cooling works. 

INTRODUCTION 

The surface resistance of the superconducting (SC) radio 
frequency (RF) cavity consists of the temperature depen¬ 
dent part and the temperature independent part. The latter 
is called the residual resistance, R res , and limits the quality 
factor of SCRF cavity at T <C T c . 

Magnetic fluxes trapped in a process of cooling down 
of a cavity degrade R res . Thus decreasing a number of 
trapped fluxes is necessary for a reduction of R res . Recent 
studies show that cooling down conditions affect a num¬ 
ber of trapped fluxes d % 3, 4]. In particular, researchers 
in Fermilab found a fast cooling with a larger temperature 
gradient leads to a better expulsion of fluxes and thus yields 
a lower residual resistance R res El SI- They achieved an 
ultra high Qo (— 2 x 10 11 ) by the fast cooling method fU. 

While many experimental studies on the flux trapping 
have been conducted, not much theoretical progress fol¬ 
lowed on it. In the present paper, we theoretically study 
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Figure 1: SC material of a cavity wall under cooling down. 
The gray and blue region represent the domains with T > 
T c and T < T c , respectively. The origin of the x-axis is 
located at the interface of these two regions, which we call 
St c • The SC material is cooled down from the right to the 
left. The ambient magnetic field B a is parallel to the z-x 
plane, and 0 is the angle between B a and the z- axis. 

the flux trapping that occurs in the process of cooling down 
by focusing on the dynamics in the vicinity of the phase 
transition fronts. We do not consider effects of the ther¬ 
mal current. We show a number of trapped fluxes and thus 
R res are proportional to the ambient magnetic field and the 
inverse of the temperature gradient. The present model re¬ 
veals why and how the residual resistance depends on the 
temperature gradient, and naturally explains how the fast 
cooling works. 

MODEL 

Let us consider the SC material shown in FigureQ] which 
represents a part of a cavity wall. The gray and blue regions 
represent the domains with T > T c and T < T c , respec¬ 
tively. The origin of the x-axis is located at the interface of 
these two regions, which we call St c in the following. The 
SC material is cooled down from the right to the left: VT is 
parallel to the negative direction of the x-axis; dT/dx < 0 
and dT/dy = dT/dz = 0. The ambient magnetic field B a 
is parallel to the z-x plane, and 0 is the angle between B a 
and the z- axis. Its magnitude |B a | is given by B a . 










Figure 2: Schematic view of the vicinity of the bound¬ 
ary St c - There exist three domains: the normal conduct¬ 
ing domain ^Normal (x < x C 2), the vortex state domain 
TVortex ( x c 2 < x < x c ±), and the Meissner state domain 
^Meissner (x > x c \). As the material is cooled down, the 
phase transition fronts sweep from the right to the left. Vor¬ 
tices contained in XVortex are trapped by pinning centers in 
this process. 


Let us look at the vicinity of St c • Since T ~ T c in the 
vicinity of St c , the lower critical field B c i(T(x)) and the 
upper critical field B C 2 (T(x)) are strongly suppressed and 
can be smaller than the ambient magnetic field B a (typi¬ 
cally ^ 0.1-10 fjb T). Then we see that there exist two types 
of the phase transition fronts in the vicinity of St c : x = x C 2 
at which B c 2 = B a and x = x c \ at which B c \ = B a . Thus 
there exist three domains as shown in Fig. [2 


^Normal • X < X c 2 
< X^Yortex • X c 2 ^ ^ X c \ 


X^Meissner • X X c \ 


( 1 ) 


where X) No rmai, ^Vortex, and X^Meissner represent the nor¬ 
mal conducting domain, the vortex state domain, and the 
Meissner state domain, respectively. As the SC material 
is cooled down, the domain XVortex with the thickness 
Sx = x c \ — x C 2 sweeps from the right to the left and fi¬ 
nally passes through the entire region of the SC material 
(see Fig. 13. In this process, vortices contained in XVortex 
are trapped by pinning centers and contribute to R res . 

EVALUATIONS OF THE NUMBER OF 
TRAPPED FLUXES AND -Rres 

The goal of this section is to evaluate a number of 
trapped fluxes and R res based on the model introduced in 
the last section. 


Review of the Temperature Dependences of Rele¬ 
vant Parameters 

Since we are interested in the vicinity of the phase transi¬ 
tion fronts, where T ~ T c , we work in the framework of the 
Ginzburg-Landau (GL) theory. Let us summarize the tem¬ 
perature dependences of relevant parameters in GL theory. 
In the following, we use the normalized temperature, 



According to GL theory, the coherence length is given 
by 


C(r) = r(i-i) _i (3) 


with = y 7 ^ 2 /|a 0 |, and the penetration depth is given 

_ 1 

by A(T) = A*(l — t'j 2 with A* = ^/3/8/ioe 2 7|cto|, 

where cto, /?, and 7 are the constants derived from the mi¬ 
croscopic theory Q. Thus n = A(T)/£(T) = A*/£*, is 
a temperature independent parameter in the framework of 
GL theory. The upper critical field is given by 


B c2 (T) = 


00 

2tt£(T) 2 


00 

2 tt £* 2 


( 1 -*), 


(4) 


and the thermodynamic critical field is given by B C (T ) = 
B c2 {T)/{V2k) = (0 0 /2v / 2tt<* 2 )(1 - t) , where 0 O = 
2.07 x 10 “ 15 Wb is the flux quantum. When n ^ 1, the 
lower critical field is given by the compact expression, 


B cl (T) = 


In k, + a 
\[2 k 


Bc(T) = 


0 o (In ft + a) 

we 2 1 j ’ 


(5) 


where a ~ 0.5. In the following, we assume n 1 and 
use Eq. © to evaluate B c \ . This assumption allows us to 
analytically grasp the physics behind the flux trapping in a 
cooling down process. 

The parameter = \j'yh 2 /\ao\, which appears in 
Eqs. (O-©, is related to the coherence length in the mi¬ 
croscopic theory, £bcs = 0). Here we give 

in two special cases as examples: = 0.74£bcs (clean 

limit) and = 0.85 V£bcs^ (dirty limit), where t is the 
electron mean free path. 


Translation into the Position Dependences 

The above parameters as functions of the temperature, 
T, can be translated into functions of the position, x, by 
using the temperature distribution, T(x). Since the origin 
of the x-axis is located at the boundary St c , where T = T c , 
the temperature distribution T{x) can be written as 


dT 

T(x)=T c +—x = T c ( i- 



( 6 ) 


Vo, /A and 7 are given by ao = — ^(0), (3 = 
7C(3)z/(0)/87t 2 kgT^, 7 = 7C(3)^(0)r;|,/487r 2 A;^T 2 for clean SCs, 
and 7 = 7riy(0)D/8hkBT c for dirty SCs, respectively, where z/(0) is the 
density of state at the Fermi energy, D is a diffusion constant, and vp is 
the Fermi velocity. 
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Figure 4: On our coordinate system moving with St c , 
the phase transition fronts, x = x C 2 and x = x c \, are at 
rest, but pinning centers move. This can be described as 
“beam-target collision” events, where beams and a target 
correspond to pinning centers and the vortex state domain, 
respectively. 


Figure 3: Examples of the temperature and the criti¬ 
cal field distributions as functions of x , where \dt/dx\ = 
T~ 1 \dT/dx\ = lm" 1 , k = 5, and £* = 20 nm are as¬ 
sumed. As an example, the ambient magnetic field B a = 
10 pT is also shown together. The regions B a > B c2 , 
B C 2 > B a > B c 1 , and B a < B cl correspond to XWmai, 
^Vortex, and ^Meissner, respectively (see also Fig. [3. 


By using Eq. © or t = 1 — \dt/dx\x , we find 1 — t = 
\dt/dx\x and (1 — t)~ 2 = 1/y/\dt/dx\x. Then the coher¬ 
ence length given by Eq. © can be written as a function of 


m*)) 



(7) 


The critical fields given by Eqs. 0 and ([5]) can also be 
written as functions of x: 


Bc2(T(x)) 

B cl {T(x)) 


0o dt 
27r£* 2 dx 
0o (In ft + a) 
4t m 2 £* 2 


dt 

dx 


( 8 ) 

(9) 


The positions of the phase transition fronts, x C 2 and x c i, 
are given by those at which B a equals B c2 and B c 1 , re¬ 
spectively. Substituting x = x C 2 (x = x c \) and B c2 = B a 
(B c i = T? a ) into Eq. ® (Eq. ©), we find 


Xc2 = 2tt r 2 


Ba 


00 



4^ 2 £* 2 B a 
In k + a 0o 


dt 


dx 


( 10 ) 

( 11 ) 


Fig. [3] shows t = T/T c , T > C 2 and T> c i as functions of 
x. The region B a > B c 2 or x < x C 2 corresponds to 
^Normal, B c2 > B a > B cl or x c2 < x < x cl corre¬ 
sponds to TVortex, and B a < B c 1 or x > x c i corresponds 

tO ^Meissner- 


It should be noted that all the above calculations are 
based on the GL theory, and are valid near T ~ T c or 
\dt/dx\x <C 1 [see also Eq. ©]. Thus the positions of the 
phase transition fronts, x c2 and x c \, given by Eqs. (fTOl) and 
O are also valid only when x c2 <^\dt/dx\ 1 and x c \ <C 
\dtfdx\~ 1 , respectively. These conditions are satisfied as 
long as a typical ambient magnetic field < 0(10) pT is 
assumed. 

Number of Trapped Fluxes 

On our coordinate system that moves with the bound¬ 
ary St c , the phase transition fronts are at rest, and pin¬ 
ning centers move and collide with vortices (see Fig. [4]). 
Then the flux trapping phenomenon can be described as a 
beam-target collision event with a reaction cross-section, cr, 
where pinning centers and the domain D vor tex correspond 
to beams and a target, respectively. 

We assume the pinning force can reach a distance of the 
order of the coherence length like that of a grain boundary 
or a normal conducting precipitate in an SC with k, 1. 
Note that the size of £ in the vicinity of the phase transition 
fronts, where T ~ T c , is much larger than that of T <C T c 
and typically 0(1)pm. For example, at the phase transition 
front x c2 , 


£(Zc2) 


r 


sfxti 


dt 

dx 



( 12 ) 


which yields £(x c2 ) ~ 6pm when B a = 10/iT, and 
£(x c2 ) ~ 60/im when B a = 0.1 pT. 

Pinning centers and vortices have effective radii ^ £. 
Since £ is maximum at x = x c2 and decreases as x in¬ 
creases, we assume reactions mostly occur near x = x c2 
for simplicity (see Fig. [4]). Then the reaction cross-section 
<7 is given by 


cr oc £(a ; c2 ) 2 . 


(13) 



























where £(x C 2 ) is given by Eq. (fl2l) . 

The number of vortices in the “target” V VO rtex is given 
by 

N i 5x) oc ^ALySx, (14) 

00 

where B a j 0o is a vortex density, and Sx is the thickness of 
domain D vor tex given by 


5x = X cl - X c2 = 47 tC fn(^) 


Ba 

00 


dt 

dx 


-l 


(15) 


and f n (ac) = ft 2 /(lnft + a) — 1/2. It should be noted 
that, even if some vortices in D vortex are trapped in the 
process of passing through the material, ’ does not de¬ 
crease, because vortices are supplied from the phase tran¬ 
sition front at x = x C 2 as long as the material is in the am¬ 
bient magnetic field. Here we give examples of the typical 
size of Sx. Assuming k = 5, = 20 nm and B a = 10 p T, 

Sx = 3 x 10 2 p m for \dt/dx\ = 1 m -1 (see also Fig. [3]) and 
Sx = 30/imfor \dt/dx\ = 10m _1 . As \dt/dx\ increases, 
the target thickness Sx and thus decrease. 

Now we can evaluate the number of trapped fluxes. 

($ x \ 

Introducing the total cross-section E = crN^ , the 
reaction probability is given by P = T>/AL y AL z . 
Then the number of events is given by A event = 
(PpinA.L x AL y AL z ^P = Ppi n ALxcrN^ \ where /9pi n is 
the density of pinning centers that have strong enough pin¬ 
ning forces to pin vortices against the forces due to the free 
energy gradient that prevent the vortex penetration. Since 
the number of trapped flux AT t rap is expected to be propor¬ 
tional to Advent, we obtain 


Ntrap OC PpinALxCrN^ OC B a 


dt 

dx 


-l 


(16) 


The factor \dtfdx\~ 1 comes from the fact that the total re¬ 
action cross-section E and thus the reaction probability P 
is proportional to the thickness of the vortex state domain, 

Sx oc \dtfdx\~ 1 . 


Residual Resistance 

Let us start from the well-known formula for the resid¬ 
ual resistance, i?re S 0% = Rn{Ba/B C 2 )- This formula as¬ 
sumes that 100 % of the ambient magnetic field is trapped. 

is reduced to R n when B a = B c 2 . This formula 
can be naturally generalized to R res = R n (r tra p B a / B C 2 ), 
where r tm p is the ratio of a number of trapped fluxes N trap 
to a number of total ambient fluxes, N^ ALx \ Note that 
R res is reduced to R^ % when r tra p = 1. Since the to¬ 
tal ambient fluxes is given by N^ ALx ^ = (A L x /Sx)N^ x \ 

We obtain rtrap = A^rap/A^, OC Ppin&Sx Or rtrap = 
M\dt/dx\~ x , where M oc p P m£* 2 /n(ft) is a material de¬ 
pendent parameter. Then we find 


Rn 


E> Arap^a _ r < D 

XCn, — C JJ a 


B t 


c2 


dt 

dx 


-l 


(17) 



Figure 5: R res as functions of dT/dx. The blue and red 
symbols represent experimental results read from Ref. EB, 
where a distance between the two sensors at the equator 
and the iris is assumed to be 0.1 m. The solid and dashed 
curves represent Eq. (l20t with B a = 1 pT and 19 p T, re¬ 
spectively. 


where C = ( R n /B c2 )M oc (R n /B c2 )p p \ n £,* 2 /»(«)• 

Put briefly, trapped fluxes have the normal cores, and the 
total normal conducting area increases as A tr a P increases 
(see also Ref. (5)): 

(18) 

As mentioned above, the factor \dt/dx\~ x comes from the 
thickness of the vortex state domain Sx. As \dT/dx\ in¬ 
creases, Sx decreases, E and P decreases, N trap decreases, 
and Rres decreases. 

It should be noted that even if we use the penetration 
depth instead of the coherence length as a distance that the 
pinning force can reach, the resultant functional form of 
Rres is unchanged because the penetration depth has the 
same temperature dependence (oc (1 — t)~ 2 ) as the coher¬ 
ence length. 

Range in Application of the Present Model 

When | dT/dx | is so large that the thickness of the vortex 
state domain Sx is smaller than £( 2 ^ 2 ), the present model 
ceases to be valid. Such \dT/dx\ is obtained by solving 
£{x C 2 ) = Sx and is given by 

= 4 7 r^* 2 /4«) iff) ■ ( 19 ) 

On the other hand, when | dT/dx | is too small, the unitarity 
is broken. Thus a cutoff exists at small \dT/dx\, where 

100 % of the ambient magnetic field is trapped: i? res = 
# 100 % 

-Ttfes 

COMPARISON WITH EXPERIMENTS 

We derived R res given by Eq. (IT71) in the last section. 
Now, we compare Eq. (TTtT) with experiments by Roma¬ 
nenko et al. a. 


dt 

dx 


Rres A^trap ^ B a 


dt 

dx 





















Let us extract the constant C from the experimental re¬ 
sults 0| . Assuming the distance between the temperature 
sensors at the equator and the iris is given by ^ 0.1 m, the 
temperature differences between the equator and iris given 
in Ref. Q) can be translated into the temperature gradients. 
Then C can be evaluated by C = \dt/dx\R res /B a , and we 
find 


R 


res 


/ 3.08 nfA dt_ 
\ m • fi T ) a dx 


( 20 ) 


Fig. [5] shows R res as functions of \dT/dx\ (= 
T c \dt/dx\). The blue and red symbols represent experi¬ 
mental results read from Ref. 0, where the temperature 
differences are translated into the temperature gradients 
by assuming the distance between the temperature sensors 
at the equator and iris is given by ^ 0.1m. The solid 
curve represents Eq. (f20l) with B a = 1 /iT (10 mG), which 
agrees well with the experimental results. The dashed 
curve represents Eq. (l20l) with B a = 19 /iT (190 mG). 
We see that, even under such a strong ambient magnetic 
field, R res ~ 5 nCl can be achieved by a large \dT/dx\ ~ 
100 K/m. This result is also consistent with Ref. 0, where 
they achieved Qo > 5 x 10 10 with the cavity cooled down 
under B a = 19 /iT (190 mG). 

Other examples and more detailed discussions will be 
presented elsewhere 0. 


SUMMARY 

In the present paper, we studied the flux trapping that 
occurs in the process of cooling down of a superconductor 
by focusing on the dynamics in the vicinity of the phase 
transition fronts. 

• We considered the simple model shown in Fig. [I] The 
gray and blue region represent the domains with T > 
T c and T < T c , respectively. The origin of the x-axis 
is located at the interface of these two regions, which 
we call St c • The SC material is cooled down from the 
right to the left. 

• In the vicinity of St c , where T ~ T c , B c \ and B c2 
are strongly suppressed and can be smaller than the 
ambient magnetic field B a ^ 0.1-10/iT. There ex¬ 
ist two phase transition fronts in the vicinity of St c : 
x = x C 2 at which B c 2 = B a and x = x c \ at which 
B c i = B a . Then there exist three domains: the nor¬ 
mal conducting domain (x < x C 2), the vortex state 
domain ( x C 2 < x < x c \), and the Meissner state do¬ 
main (x > x c i). See Fig. [2] and [3] 

• As the material is cooled down, the phase transition 
fronts together with the vortex state domain sweep the 
material. Vortices contained in the vortex state do¬ 
main are trapped by pinning centers in this process 
(see Fig.0. 

• A number of trapped fluxes iVt rap were evaluated by 
using the analogy with the beam-target collision event. 


Beams and a target correspond to pinning centers and 
vortex state domain, respectively (see Fig. 0 ). 

• The number of trapped fluxes 7Vt rap is given by 
Eq. (IT6l) . which is proportional to B a and \dTfdx\~ 1 . 

• Trapped fluxes have normal cores, and the total nor¬ 
mal conducting area increases as N trap increases. 
Thus Rres is proportional to N trap and is given 
by Eq. £0, which is also proportional to B a and 

\dTfdx\~ 1 . 

• The factor \dTfdx\~ 1 in N trap and R res comes from 
the fact that the thickness of the vortex domain is pro¬ 
portional to \dT/dx\~ lm . As \dT/dx\ increases, the 
thickness of the vortex state domain decreases, the to¬ 
tal reaction cross-section and the reaction probability 
decrease, then a number of trapped fluxes N trap de¬ 
creases, and the residual resistance, R res decreases. 

• The residual resistance, R res , given by Eq. (flTt was 
compared with experimental results (see Fig. 0 . The 
blue and red symbols represent experimental results 
read from Ref. 0 , where the temperature differences 
in the original paper are translated into the temper¬ 
ature gradients. The solid curve represents Eq. (l20l) 
with B a = 1 /iT, which agrees well with the experi¬ 
mental results. Note that our formula for R res contains 
only one free parameter. 

• The dashed curve in Fig. [5] represents Eq. (l20t with 
B a = 19 /iT (190 mG). We see that, even under such 
a strong ambient magnetic field, R res ^ 5ncan be 
achieved by a large \dT/dx\ ~ 100 K/m. This result 
is also consistent with Ref. E), where they achieved 
Qo > 5 x 10 10 with the cavity cooled down under 
B a = 19 /iT (190 mG). 

The present model revealed why and how R res depends 
on the temperature gradient, and naturally explained how 
the fast cooling works. Other examples and more detailed 
discussions will be presented elsewhere 0. 
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